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We consider a dilute granular gas of hard spheres colliding inelastically with coefficients of normal 
and tangential restitution a and j3, respectively. The basic quantities characterizing the distribution 
function /(v, bj) of linear (v) and angular (u>) velocities are the second-degree moments defining 
the translational (T tr ) and rotational (T Iot ) temperatures. The deviation of / from the Maxwellian 
distribution parameterized by T tr and T rot can be measured by the cumulants associated with the 
fourth-degree velocity moments. The main objective of this paper is the evaluation of the collisional 
rates of change of these second- and fourth-degree moments by means of a Sonine approximation. 
The results are subsequently applied to the computation of the temperature ratio j 11 " '/! 1 ' 1 anc i the 
cumulants of two paradigmatic states: the homogeneous cooling state and the homogeneous steady 
state driven by a white-noise stochastic thermostat. It is found in both cases that the Maxwellian 
approximation for the temperature ratio does not deviate much from the Sonine prediction. On the 
other hand, non-Maxwellian properties measured by the cumulants cannot be ignored, especially 
in the homogeneous cooling state for medium and small roughness. In that state, moreover, the 
cumulant directly related to the translational velocity differs in the quasi-smooth limit /3 — > — 1 from 
that of pure smooth spheres (f3 — —1). This singular behavior is directly related to the unsteady 
character of the homogeneous cooling state and thus it is absent in the stochastic thermostat case. 



I. INTRODUCTION 

Among the many topics in the kinetic theory of gases 
uncovered by Carlo Cercignani during his long and fruit- 
ful scientific career it is mandatory to mention the kinetic 
theory of inelastic particles, a field he substantially con- 
tributed to during the last decade of his With this 
paper we wish to pay a modest tribute to Carlo Cercig- 
nani's accomplishments in this field. 

The most frequently used physical model of a granular 
fluid consists of a system of many inelastic and smooth 
hard spheres with a constant coefficient of normal resti- 
tution a.— On the other hand, the macroscopic nature of 
the grains makes the influence of friction when two par- 
ticles collide practically unavoidable*^— From a more 
fundamental point of view, the existence of collisional 
friction is important to unveil the inherent breakdown of 
energy equipartition in granular fluids, even in homoge- 
neous and isotropic states. 

The simplest model accounting for friction during colli- 
sions assumes, apart from a constant coefficient of normal 
restitution a, a constant coefficient of tangential restitu- 
tion /3i 12 ' 13 While a is a positive quantity smaller than 
or equal to 1 (the value a = 1 corresponding to elastic 
spheres), the parameter f3 lies in the range between —1 
(perfectly smooth spheres) to 1 (perfectly rough spheres). 
The total kinetic energy is not conserved in a collision, 
unless a = 1 and (3 = ±1. As a consequence, many of 
the papers in the literature assume that the spheres are 
nearly smooth and nearly elasticJ^— 

The theoretical study of a granular gas is usually 
undertaken by employing tools already developed in 



nonequilibrium statistical mechanics and kinetic theory 
of normal gases. In particular, one can introduce the 
one-body distribution function /(r, v, w;i), where v and 
u> are the velocity of the center of mass and the angular 
velocity, respectively, of a particle. From the second- 
degree velocity moments of the distribution function it 
is straightforward to define (granular) translational and 
rotational temperatures, T tT and T rot (see Sec. [TTJ) . The 
rates of change of these two quantities produced by col- 
lisions define the energy production rates £ tr and £ rot as 

1 /<9T tr \ „ ot 1 fdT lot " 



T tr V dt 



^ rot V dt / con 

(1.1) 

The collisional energy production rates £ tr and £ rot do 
not have a definite sign. They can be decomposed into 
two classes of terms^ equipartition rates and cooling 
rates (see Fig. []}. The equipartition terms, which exist 
even when energy is conserved by collisions (a = 1 and 
j3 = ±1), tend to make temperatures equah 35 ' 37 . There- 
fore, they can be positive or negative depending essen- 
tially on the sign of the temperature difference r r tI — r r TOt . 
On the other hand, the genuine cooling terms reflect the 
collisional energy dissipation and thus they are positive 
if a < 1 and/or \(3\ < 1, vanishing otherwise. Only the 
cooling terms in £ tr and £ rot contribute to the net cool- 
ing rate C = (£ tr T tr + <f ° t T rot )/(T tr + T rot ). Both £ tr 
and £ rot are functional of / and therefore they formally 
depend on all the moments of /, not just on T tr and T Iot . 

In an extensive paper^ Goldshtein and Shapiro un- 
dertook the task of evaluating the collisional energy pro- 
duction rates £ tr and £ rot by using a two-temperature 
Maxwellian approximation for the distribution function, 
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FIG. 1: (Color online) Scheme on the two classes of contri- 
butions (equipartition rates and cooling rates) to the energy 
production rates £ tr and £ rot characterizing the effect of col- 
lisions on T tr and T Iot , respectively. The terms represented 
by dotted arrows are absent in the case of perfectly smooth 
spheres (j3 = — 1). 
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Here, n and u are the number density and the flow ve- 
locity, respectively, of the gas and m and / are the mass 
and the moment of inertia, respectively, of a particle. 
The mean angular velocity has been assumed to vanish.— 
The final results for the energy production rates in the 
Maxwellian approximation ar o 29 ' 34 
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v, (1-4) 
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is the rotational/translational temperature ratio, 
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is the dimensionless moment of inertia {a being the di- 
ameter of a particle) , and 
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v = —a 2 n^yirT tT /m 



(1.7) 



is an effective collision frequency. The expressions within 
the Maxwellian approximation but with a non-zero mean 
angular velocity can be found in Ref. [36]. Furthermore, 
the more general expressions for mixtures were derived 
in Ref. M- 

Equations (|1.3[) and (|1.4[) have been applied to the so- 
called homogeneous cooling state (HCS)i 22 i 26 From the 



condition lim^oo 9(t) = const one gets £ tr = £ rot , yield- 
ing a quadratic equation for the asymptotic temperature 
ratio 8 whose physical solution is 



9 = \/l + C 2 + C 



(1.8) 



with 



C = 



2/e(l + j8) 
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or 
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(1.9) 

The time evolution of the ratio T rot /T tr toward the HCS 
asymptotic value (|1.8p has been widely analyzed, both 
theoretically and by means of molecular dynamics, by 
Luding, Zippclius, and co-workers 

An even simpler application of Eqs. (|1.3[) and (|1.4[) 
corresponds to the case of a homogeneous and isotropic 
granular gas kept in a nonequilibrium steady state by a 
white-noise thermostat <22r— We will refer to this situa- 
tion as the white- noise state (WNS). The steady-state 
condition £ rot = simply yields 
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(1.10) 



Despite the crudencss of the Maxwellian approxima- 
tion given by Eq. (|1.2p . Eq. (|1.9p does a very good 
job when compared with computer simulations for the 
HCS The same is expected to hold for Eq. ([TTUj) in 
the WNS case. On the other hand, the production rates 
£ tr and £ rot , being nonlinear functionals of /, can be ex- 
pected to be influenced by non-Maxwellian features of 
/, thus deviating (even if only slightly) from Eqs. (|1.3|) 
and (|1.4[) . The basic non-Maxwellian features of a ve- 
locity distribution function /(v,u>) are the existence of 
non-zero cumulants. The most physically interesting cu- 
mulants are a2o, flu, and 0,02, defined as 
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(1 + au), (1-12) 
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(1 + 002). (1-13) 



Here the angular brackets denote average values defined 

as 

(V(v,«)> = - fdv [ dwV(v,a»)/(v,w). (1.14) 



The objectives of this paper are: (a) to evaluate the 
second-degree collisional moments £ tr and £ rot in a Soninc 
approximation that includes the cumulants a 2 o, an, and 
0-02; (b) to evaluate the three fourth-degree collisional 
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moments related to the moments denned by Eqs. 
(|1.13|) in the same Sonine approximation; and (c) to ap- 
ply the results to both the HCS and the WNS in order to 
"refine" Eqs. (|1.8[) and (|1.10[) , and estimate 020, an, and 
ao2 in those states. The method will be similar to that 
already worked out in the case of smooth spheres.— ~— 

This paper is organized as follows. The collision rules 
and the Boltzmann equation for a gas of inelastic and 
rough hard spheres are presented in Sec. |TTJ The Sonine 
approximation is constructed in Sec. IIIH where the de- 
rived expressions for the collisional moments are written 
down. Sections IIVI and [V] deal with the application of 
the results to the HCS and the WNS, respectively. The 
paper ends with a brief discussion in Sec. IVII 



II. 



COLLISION RULES AND BOLTZMANN 
EQUATION 

A. Collision rules 



Let us consider a granular gas made of inelastic rough 
hard spheres of mass m, diameter a, and moment of in- 
ertia /. In this section we first derive the rules for a 
binary collision between two spheres with precollisional 
center-of-mass velocities (vi,V2) and angular velocities 
(wi,w 2 ). 

Let us denote by V12 = vi — V2 the precollisional rela- 
tive velocity of the center of mass of both spheres and by 
a = (r 2 — i"i)/|r 2 — i"i| the unit vector pointing from the 
center of sphere 1 to the center of sphere 2. The precol- 
lisional velocities of the points of the spheres which are 
in contact during the collision are 

Vi =vi-^5xu,. V2 = v 2 l^xu 2 , (2.1) 

the corresponding relative velocity being 

V12 = V12 - a x S12, S i2 = ~(u>: + w 2 ). (2.2) 



Conservation of linear and angular momenta yield: 
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v[ + v 2 = vi + v 2 , (2.3) 



— m -^ cr x v i = It*>i — m 7£ <T x v i; (2.4a) 



Iu)' 2 + x v 2 = + m l^ a x v 2, (2-4b) 

where the primes denote postcollisional values. Equa- 
tions (|2~3)) and (j2~4)l imply that 

vi=Vi-A 12 , v 2 = v 2 + Ai2, (2.5) 



, ma ^ A , ma ^ 

= wi - —a x A12, lv 2 =uj 2 -— crx A i2 , (2.6) 



where mAi 2 is the impulse exerted by particle 1 on par- 
ticle 2. Therefore, 

2 

v 12 = V12-2A12, Vi 2 = Vi 2 -2Ai 2 H — ffx (5- x A X2 ) , 

K 

(2.7) 

where the dimensionless moment of inertia k is defined 
by Eq. (|1.6[) . It varies from zero to a maximum value 
of |, the former corresponding to a concentration of the 
mass at the center of the sphere, while the latter value 
corresponds to a concentration of the mass on the surface 
of the sphere. The value k = | refers to a uniform mass 
distribution. 

To close the collision rules, we need to express A 12 in 
terms of the precollisional velocities and the unit vector 
<t. To that end, let us relate the normal (i.e., parallel to 
tr) and tangential (i.e., orthogonal to <x) components of 
the relative velocities V12 and V' 12 by 

a ■ V' 12 = -aa ■ V12, a- x V' 12 = -fta x Vi 2 . (2.8) 

Here, as said in Sec. HI a and 8 are the coefficients of nor- 
mal and tangential restitution, respectively. The former 
coefficient ranges from a — (perfectly inelastic parti- 
cles) to a = 1 (perfectly clastic particles) , while the latter 
runs from /3 = — 1 (perfectly smooth particles) to = 1 
(perfectly rough particles). A more realistic model con- 
sists of assuming that (3 is a function of the angle between 
V12 and <x^ thus accounting for Coulomb friction. In 
this paper, however, we will assume a constant (3. 

Inserting the second equality of Eq. (|2.7p into Eq. (|2.8[) 
one gets 

<r • A12 = acr ■ V12, <r x A12 = PS x V12, (2.9) 
where the following abbreviations were introduced: 

_ 1 + a ~ k 1 + 8 

a=——, B=- --. 2.10 

2 1+k 2 v ; 



Therefore, 

A12 = 5(vi 2 • er)er + 8 [V12 - (V12 • ct)(t -ax S i2 ] . 

(2.11) 

Equations (|2.5[) . (|2.6[) . and (|2 . 1 L |) express the postcol- 
lisional velocities (v^wi, v 2 ,w 2 ) in terms of the precol- 
lisional velocities (vi, u}\, V2, w 2 ) and the unit vector <r. 
In the special case of perfectly smooth spheres (f3 = — 1 
or, equivalently, /3 = 0) one has 5 x A 12 = 0, so that 
uj[ = u>i and U3 2 = u 2 . 

The collisional change of the total (translational plus 
rotational) kinetic energy is 

E' l2 -E 12 = ~ (1 - a 2 ) (a ■ v 12 ) 2 

(1 - /3 2 ) [via 



4 1 + K 



-SxS 12 -(v 12 -?)5f, (2.12) 



where 



E\2 



2 1 2 , 2 



(2.13) 



4 



The right-hand side of Eq. (|2.12[) is a negative definite 
quantity. Thus, energy is conserved only if the particles 
are elastic (a = 1) and either perfectly smooth (/3 = —1) 
or perfectly rough (f3 = 1). Otherwise, E' 12 < E\ 2 and 
kinetic energy is dissipated upon collisions. 

Equations (|2.5[) . (|2.6p . and (|2. 1 1|) give the direct col- 
lisions! rules. For a restituting encounter the pre- and 
postcollisional velocities are denoted by (v", uj'{, v 2 ' , u> 2 ) 
and (vi, oj\, \ 2 , u> 2 ), respectively, and the collision vector 
by a" = —a. It is easy to verify that the relationship 
V12 • a" = —av'{ 2 ■ = — V12 • a holds. Analogously, 
a" x V12 = —fta" x V12 = — & x Vi2- As a consequence, 
the restituting collision rules are 



v" = V! - A l; 



= v 2 + A12, 



(2.14) 



u>x-— crxAi2, u 2 = u 2 -— ax A12, (2.15) 
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where 

A12 



(v 12 • a)a -ax S 12 ] 



a _ /3 

-(Vi2 ' <T)cr + - V12 

a p 

(2.16) 

The modulus of the Jacobian of the transformation 



between pre- and postcollisional velocities is 



d(vi,u>i, V 2 ,U> 2 ) 



d(\i, wi, V2,u> 2 ) 



a/3 2 



(2.17) 

Furthermore, the relationship between volume elements 
in velocity space reads 



\v'{ 2 -a"\dv'{du>1dv%du> 2 ' 



|Vi2 ■ <r| 

a 2 /3 2 



dvxd<jJidv2ckiJ2. 



(2.18) 



B. Boltzmann equation 



If the granular gas is dilute enough the velocity 
distribution function /(r,v,u;;£) obeys the Boltzmann 



,22.44 



cquatior 

a t / + v-v/ = j[v,w|/], 

where the collision operator is 



(2.19) 



J[vi,wi|/] = a 2 J dv 2 J du> 2 J da (v 12 • ?) 



a 2 /3 : 



r/i /2 _ /1/2 



(2.20) 



Here the subscript + in the integral over a means the 
constraint V12 • a > and we have employed the short- 
hand notation /{' = /(v'/, u>") and so on. 

Given an arbitrary function ^>(v, a;), its average value 
is defined by Eq. (|1.14j) . The associated collisional rate 



of change is (ip) 1 J'[ip\f], where the collisional quantity 
J[^}{v,Lj)\f\ is defined by 



Jim 



J d\\ J dux -0(vi,u;i)J[vi,u;i|/] 



dvi J du>i J dv 2 J duj 2 J da 
x (vis • a)fxf 2 (ip[ +ip' 2 -tpi- ip 2 ) , 



(2.21) 

where in the last step we have carried out a standard 
change of variables. 

Here we are especially concerned with the partial tem- 
peratures associated with the translational and rotational 
degrees of freedom: 



T 



-<(v-u) 2 ), r- 



(2.22) 



where u = (v) is the flow velocity. The corresponding 
energy production rates are defined as 



-j[( V -uy\f], e 



3nT 



^Jl"V}- 
(2.23) 



The total temperature and its corresponding cooling rate 
are 



T = 



rptr _|_ ^rot 



(2.24) 



(2.25) 



It is 

of T rot , 



worthwhile remarking that, instead 
we could have alternatively adopted 
(J/3)((w- (w)) 2 ) = T Iot (l-X), with 
X = Kmcr 2 (u;) 2 /12T rot , as the definition of the ro- 
tational temperature. However, a disadvantage of this 
alternative choice is that, in contrast to the cooling 
rate £ defined by Eq. (|2.25[) . the alternative "cooling" 
rate £ associated with the alternative total temperature 
T = i(T tr + T rot ) = T — \T mt X is not positive definite 
and in fact becomes negative in the perfectly elastic and 
rough case (a = 1, /3 = l)i^£ 

Making use of the collision rules given by Eqs. (|2.5I) . 
(|2.6[) . and (|2.11[) , and after performing the integration 
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over a. one gets 

tr 50TI/ 



96(T tr /m) 3 / 2 
—pr{{3vi 2 S 2 2 



(5(1 -a) +0(1-0)) 



2 (V12 • s 12 y 



(2.26) 



96(T tr /m) 3 / 2 

~Vl 2 ( V 12 • S12)' 



P, 



J \ 2 i 



(2.27) 
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384(T t V?7i) 3 / 2 1 
k 1-/3 2 



(1-« 2 )K 2 » 



1 + K 2 
-V12 ( v 12 ' Sl 2 ) 2 



2(K 2 }) + ((touSl 



(2.28) 



In these equations, 9 is the temperature ratio defined by 
Eq. (|1.5j) . 1/ is the collision frequency defined by Eq. (|1.7jl . 
and 



<KV12,S 12 ))> 



1 



Xf/;(vi2,Si2)/(vi,u;i)/(v2,u;2) 

(2.29) 



are two-body averages. Use has been made of Eq. (|2.10j) 
upon obtaining Eq. (f2T2"g| from Eqs. ([23d) and ([2~27]) . 
It is worthwhile emphasizing that Eqs. (|2.26[) f|2.28[) arc 
exact in the framework of the Boltzmann equation. 



III. SONINE APPROXIMATION FOR SECOND- 
AND FOURTH-DEGREE COLLISIONAL 
MOMENTS 



Equations (|2.26[) and f|2 .27[) express the translational 
and rotational energy production rates as functionals of / 
through two independent two-body averages of the form 
given by Eq. (|2.29[) . If / is replaced by the Maxwellian 
approximation Eq. (|1.2[) one gets Eqs. (|1.3[) and (| 1 .4|) . As 
said in Sec. U we want to go beyond such a Maxwellian 
approximation. 

To proceed, it is convenient to introduce the dimen- 
sionlcss velocities 



u 



y/2T tr /m ^/2T" ot /I 
and the dimensionless distribution function 

1 / ^ptrrpYot \ 3/2 



(c,w) 



ml 



/(v,w). 



(3.1) 



(3.2) 



In terms of the reduced translational and rotational ve- 
locities, the collision rules given by Eqs. (|2.5p . (|2.6p . and 
(|2~TT|) become 



Wj = Wi — 



c 2 + AI 2 , 



(3.3) 



w 2 



<k6 



a x A* 12 , 



(3.4) 



A^a = a (C12 ■ <r) <t + j8 



C12 - (C12 • tr) tr 



er x (wi + w 2 ) 



(3.5) 



Let us now specialize to isotropic states. The lat- 
ter condition implies that the scalar function <f>{c, w) 
is invariant under orthogonal transformations, including 
those with determinant equal to +1 (rotations) or — 1 
(reflections). This means that <^>(c,w) is actually a func- 
tion of the three scalar quantities c 2 = c • c, w 2 = w • w, 
and (c • w) 2 . We do not need to assume that the state 
is either homogeneous or stationary. Here we focus on 
the following second- and fourth-degree moments: (c 2 ), 
(w 2 ), (c 4 ), (c 2 w 2 ), and (w 4 ). By construction, 



(c 2 ) = (w 2 ) 



In the Maxwellian approximation (jl.2l) . i.e., 
0(c,w) -> 0Af(c,w) = 7r _3 e _c2_t " 2 



one has 



15 

T' 



(c 2 w 2 ) -> 



15 



(3.6) 



(3.7) 



(3.8) 



In general, however, cj> ^ 4>m and the above equalities are 
not verified. This can be characterized by the cumulants 



a-20 



15 



1. 



«H = g^W 2 ) - 1, 



«02 = -^( W ) - L 



(3.9) 



(3.10) 



(3.11) 



Note that Eqs. PT^ )) — (|5TTT]| arc equivalent to Eqs. (fTTTTj) 

una). 

Let us define the collisional moments [i pq (with p, q = 
even) as 



J dc J dwcPw q J*[c,w\(f)], (3.12) 



where the dimensionless collision operator J* is defined 
similarly to Eq. (|2.20[) . except that one must formally 
take a — 1 and the collision rules arc given by Eqs. 
p.3p - (|3.5[) . The energy production rates £ tr and £ rot are 
directly related to the collisional moments /U20 and /io2 
by 



e = 



5v 



=M2o, r 



5// 



12V2tt 12V2tt 
Analogously, the total cooling rate is 
hv ^20 + ^02 



M02- 



nVz^ 1 



(3.13) 



(3.14) 



The primary objective in this section is to get esti- 
mates of the second-degree collisional moments )i2o and 
/102, and of the fourth-degree collisional moments (14$, 
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/122, and ^04 i n terms of the temperature ratio 9 and the 
cumulants a 2 o, an, and 002- To that end, we first express 
the distribution function <fi by the first few terms in its 
Sonine expansion, 

0(c,w) « ^/(c,w) [l + a 20 5 ( 1 2) (c 2 ) + a 02 M 2) (^) 

2 2 

+a 11 5 ( 1 1) (c 2 )5 ( 1 1) ( W 2 )" . 

2 2 

(3.15) 



The Sonine polynomials in Eq. (|3.15p arc 

M 1} (a?) = | - a?, M 2) (a?) = £ (15 - 20.t + 4x 2 ) . 

2 2 2 8 

(3.16) 

In principle, apart from the moments (c 4 ), (c 2 w 2 ), and 
(ic 4 ), the other independent fourth-degree moment ((c • 
w) 2 ) should be represented in the truncated expansion 
(13.151) . However, for simplicity, it is assumed here that 



((c-w) 2 



This implies that the study of the orientational corre- 
lation between c and w is not addressed in this paper. 
From that point of view, our approach is complemen- 
tary to that of Refs. (HH where it was assumed that 
«20 = an = ao2 = but ((c ■ w) 2 )/(c 2 u) 2 ) ^ ±. 

The second step consists of inserting the approxima- 
tion defined by Eq. (|3.15[) into Eq. (|3. 12|) , and neglecting 
terms nonlinear in 020, an, and 002- After some algebra 
one gets the following expressions for the second-degree 
collisional moments: 



^20 



4V2tt 



5(1 -5) +0(1-0) 1 + 



3a 



2(1 



16 



. 6 t (i_£?° + ^I 
k \ 16 4 



(3.18) 



M02 



/3 



I 1 - - I (l - — + — " 
I k J \ 16 4 j 



3a2o\ 
16 ) 



(3.19) 



Thus, Eq. (|3~T4|) gives 
5v 



12(1 + 6*) 

1-/3 2 



-0- 



1 



1 



1 - 



- a 
16 



K- 



1 + K 



3o2o\ 

16 J 



an 
4 



(3.20) 



Equations (|3.18[) - (|3.20[) can also be obtained from Eqs. 
(j236]) - (|2~S8f by taking into account that 
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3a 2 
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(3.17) ((3v 12 S 2 2 - v^ 1 (v 12 • Si 5 



x 1 



)• 


(3.21) 


/2T tr 


^3/2 


V m 




a 2 o 


Oll\ 


16 
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(3.22) 



in the Sonine approximation p,15p . This explains why 
the cumulant 0-02, being related to (w 4 ), does not inter- 
vene in Eqs. (f37T5|) - (|3T20|) . 

Of course, Eqs. (|3~1"8)) and (|3T9)l reduce to Eqs. ([TT5]) 
and (|1.4j) . respectively, by setting 020 = an = 0. More- 
over, Eq. p,18[) is consistent with van Noije and Ernst's 
derivation 41 - for the smooth case ((3 = —1). 

The evaluation of the fourth-degree collisional mo- 
ments /^4o, (J-22, and /io4 is much more involved. After 
carefully performing the calculations following several in- 
dependent routes to check the results, we have found 
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M40 = 16V27T <^ a A (2 - 5) + /3 d (2 - /?) - 5/3(1 - 5 - (3 + 5/3) + —(5 + 0) - —(a 2 + f3 2 ) 
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To the best of our knowledge, the collisional moments 
/Mo, M22, and /xo4 have not been evaluated before, even 
in the Maxwcllian approximation (a 2 o = an = a 2 = 0). 
The only exception is van Noijc and Ernst's evaluation 
of {I4Q in the smooth case; 41 - to which Eq. (|3.23[) reduces 
by setting /? = — 1. As an additional simple consistency 
test, we get 1122 = §/^2o m the special case of smooth 
spheres {fi = — 1) with an = 0. 

Equations (|3~T8| . f339")) . and (|3T23|) f|3T25|) are the 
main results of this paper. In the next two sections they 
are applied to the HCS and the WNS. 



IV. APPLICATION TO THE HOMOGENEOUS 
COOLING STATE 

In the homogeneous free cooling state (HCS) the Boltz- 
mann equation (|2.19|) becomes 



a t /(v,w;t) = J[v,w|/]. 



(4.1) 



As a consequence, the only mechanisms responsible for 
changes in the partial and total temperatures arc colli- 
sions. More specifically, 



-£ tr T tr , <9 t T r 



d t T = -(T. 



^■rot^rot 



(4.2) 



(4.3) 



The evolution equation for the temperature ratio 9 = 

T roy T tr if , g t Q = _^rot _ £tr)0_ 

Carrying out the change to dimcnsionlcss variables de- 
fined by Eqs. ([33]) and (j3~2")l . Eq. @3J can be rewritten 
as 

ds4>+ ^flk ' (C0) + " (W0) = J * [C ' W|0] ' (4 ' 4) 



where 9 S = (na 2 -J '2T tr / m) x 9t and use has been made 
of Eq. (|3.13j) . Taking moments in Eq. (|4~4]) we get 



-ds(<?w q ) + -(p// 20 + q^m)^) 



(4.5) 



After a certain transient period, it is expected that 
the system reaches an asymptotic regime where all the 
time dependence of / appears through one temperature 
(say T tr ) and the temperature ratio 9 remains constant. 
This implies a similarity solution of Eq. (|4.1[) of the form 
given by Eq. (|3.2j) with d s <f> = 0. Moreover, the condition 
d t 9 = implies £ tr = £ Iot = ( or, equivalently, 



M20 = M02- 

In the asymptotic regime, Eq. (I4.5[) yields 

M40 



5^2 



1 + a 2 o 



7j(M20 + M02) 



M22 
1 + an 



02 



M04 



1 + ao2 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



The objective now is to estimate the temperature ra- 
tio 9 and the cumulants 020, an, and 002 in the HCS. To 
that end, we insert the approximate expressions given 
by Eqs. (|3T8)l . ([339]) . and ([3T23 ]) -([3T25 ) into Eqs. (j4U)) - 
(|4.9[) , neglecting again terms nonlinear in 020, an, and 
ao2- Note that, for instance, Eq. (|4.7p could also be writ- 
ten as 5/i2o(l + 020) = M40>^ However, the linearization 
process gives a result different from the one obtained from 
the form (|4.7[) , as discussed in Ref . |43|. We have chosen 
the route (|4.7[) because it yields results more accurate for 
smooth spheres than the other one i 42 ' 43 
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FIG. 2: (Color online) Plot of the HCS temperature ratio 
T ro jT r vs the coefficient of normal restitution a for /3 = 
—0.9 (top panel), f3 = (middle panel), and /3 = 0.9 (bottom 
panel) . The inset in the top panel is a blow-up of the region 
0.95 < a < 1. The dashed and solid lines are the Maxwellian 
and Sonine approximations, respectively. 



From the linearized versions of Eqs. (|4.6[) — (|4.8[) one can 
obtain the three cumulants 0,20, an, and 002 as nonlinear 
functions of a, f3, and 6. Insertion into Eq. (|4.9p yields 
an eighth-degree equation for (9, whose physical solution 
is chosen as the one close to the solution (jl.8l) in the 
Maxwellian approximation. The final expressions are too 
cumbersome to be explicitly reproduced here but they are 
easy to deal with the help of a computer algebra system. 

To illustrate the dependence of 6, 020, an, and 002 on 
both a and /3, we plot those quantities as functions of a 
for three representative values of the coefficient of tan- 
gential restitution: f3 = —0.9 (small roughness), (3 = 
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FIG. 3: (Color online) Plot of the HCS cumulants 020 (top 
panel), an (middle panel), and (Z02 (bottom panel) vs the 
coefficient of normal restitution a for j3 = —0.9 (solid lines), 
ft = (dashed lines), and /3 = 0.9 (dash-dotted lines). 



(medium roughness), and (3 = 0.9 (large roughness). In 
all the cases the density of the spheres has been assumed 
to be uniform, so that k = |. The results are displayed in 
Figs. [U and [3] In Fig. [2] we observe that the Maxwellian 
approximation [cf. Eq. (|1.8jl ] does an excellent job in 
estimating the temperature ratio j llot /T tr , as confirmed 
by simulations ! 25 ' 29 This is especially true for both small 
(/3 = —0.9) and large (/3 = 0.9) roughness. While the 
Maxwellian approximation overestimates the tempera- 
ture ratio T rot /T tr for [3 = -0.9 and [3 = 0, it slightly 
underestimates this ratio for (3 = 0.9. It is interesting 
to note that, for nearly smooth spheres (/3 = —0.9), 
T rot /T tr abruptly changes from small values for nearly 
clastic spheres (a > 0.97) to very large values for in- 
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(Color online) 

\2r7-irot Irj-itr 



Plot of the HCS limit value 



FIG. 4 

lim J 9^._i(l + / g) 2 r rot /T" vs the coefficient of normal resti- 
tution a. The dashed and solid lines are the Maxwellian and 
Sonine approximations, respectively. 



clastic spheres (a < 0.95). This effect becomes more 
and more dramatic as one approaches the smooth-sphere 
limit (/3 -> -1)^ 

As for the cumulants, Fig. [3] shows some interesting 
features. In general, for large roughness (j3 = 0.9) the 
magnitudes of the three cumulants are relatively small, 
meaning that the velocity distribution function is not far 
from a Maxwellian. This agrees with the almost indis- 
tinguishability between the Maxwellian and Sonine ap- 
proximations observed in the bottom panel of Fig. [2] 
For medium roughness (fi — 0) and large inelasticity 
(a < 0.9), however, the cumulants reach relatively im- 
portant values, especially in the case of an. This trend 
is continued as roughness decreases (j3 = —0.9), except 
in the case of ao2- The latter quantity takes a high maxi- 
mum value at a ~ 0.95. This value is even higher than 1, 
thus invalidating (at a quantitative level) the lineariza- 
tion method followed to estimate it. In any case, we 
expect that the Sonine method employed in this paper 
captures the main qualitative behavior of the cumulants 
for j3 = —0.9 and a ~ 0.95. The peculiar change in the 
behavior of the cumulants when going from inelastic to 
nearly elastic spheres for small roughness (J3 — —0.9) is 
correlated to the one observed in the case of the temper- 
ature ratio. 

The smooth-sphere limit [3 — > — 1 (or, equivalently, 
/? — > 0) deserves a separate treatment. According to 
Eq. (jl.8p . one can expect in that limits the asymptotic 
behaviors 



FIG. 5: (Color online) Plot of the HCS cumulant a2o vs the 
coefficient of normal restitution a for p = — 1 (dashed line) 
and in the limit p — >■ — 1 (solid line). 



(|3.25p . and taking the limit /3 — > 0, one gets 
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(4.13) 
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16 
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3 V 16 2 8 J 



(4-14) 




(4.10) 



with di = 5(1 — a) and $2 = 1/(1 — «) in the Maxwellian 
approximation. Let us assume that a < 1. Thus, in- 
serting 9 = i} 1 Kf3- 2 into Eqs. (|3TT5j) . ([335]) . and 1(535]) - 



^04 = 20^- f 1 



a 20 Oil 
"l6 T" 



002 



(4-15) 



Since ^ 02 = £>(/?), Eqs. (|H) ]) -(prg |) imply that, in the 
smooth-sphere limit, 



M20 = 0, Hao = 0, /Li 2 2 = 0. 



(4.16) 
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Substitution of Eqs. (|4~T2|) and (|47l5|) into Eq. and 
neglecting nonlinear terms, gives an = 0. Next, from 
Eqs. (|4.11j) and (|4.14j) . together with ^ 20 = M22 = 0, we 
obtain 



a 20 = 16 



fli - 5(1 - a) 
+ 35(1-5)' 



(4.17) 



activated when the rotational temperature is sufficiently 
larger than the translational one. In other words, even if 
the dotted arrows in Fig. Q] are very weak, the great dis- 
parity between T rot and T tr triggers the flux of energy 
from the rotational toward the translational degrees of 
freedom, thus significantly modifying the cumulant a 2 o 
with respect to the case of strict smooth spheres. 



«02 = --r 



8 5(1-5) 



5 i?i +35(1-5)' 



(4.18) 



Finally, from Eqs. ([435]) . ([4"TT) . and (|4Tl8|) . together 
with the condition /140 = 0, we get a closed quadratic 
equation for §1, 

[65 - 85(1 - S)]0? - 405[9 - 125 + 45 2 (2-5)]tfi 
+ 55 2 (l-5)[59- 755 + 245 2 (2-5)] =0.(4.19) 

Figure ED shows lim^_i(l + (3) 2 T mt /T tr = [4(1 + 
k) 2 /k]i?i, where $1 is the physical solution of Eq. (|4.19D . 
as a function of a. We observe again a very good agree- 
ment between the Maxwellian and the Sonine approxi- 
mations. 

Insertion of the solution of Eq. (|4.19[) into Eqs. (|4.17[) 
and (|4. 1 8[) gives a 2 o and ao 2 , respectively. In the elastic 
limit (a —> 1) Eq. (14. 19)) yields #1 — >• 1 — 5, what implies 
fl20 — > and ao 2 — > — §. In fact, $1 ps 5(1 — 5) even in 
the inelastic case, so that ao 2 ~ — | for all a. 

The cumulant lim i g_ v _i a 2 o is plotted in Figure [5] For 
comparison, this figure also includes the curve represent- 
ing a 2 o in the pure smooth case (/3 = — 1 from the very 
beginning). In the latter case, the translational and rota- 
tional degrees of freedom are absolutely decoupled and, 
in addition, each particle keeps its initial angular velocity, 
so the arbitrary initial distribution of angular velocities 
does not change with time. This implies that only Eq. 
(|4.7[) keeps being meaningful if f3 = — 1. Inserting Eqs. 
([3~T5j) and (535)1 with f3 = into Eq. (£7) one getsi 2 -^ 



020 = -16(1 - a] 



16(1 -a) 



63 



1 - 85(1 - 5) 
- 235-85 2 (2 - 
l-2a 2 



5) 



97- 33a-2a 2 (l-a)' 



(4.20) 



We observe that a 2 o at j3 = —1 differs from lim^_j._i a 2 o. 
This singular effect of a 20 is analogous to the one ob- 
served for the ratio ((c • w) 2 ) / (c 2 w 2 ) ) 32 i 33 as well as in 
the case of the translational/translational temperature 
ratio in mixtures i 34 ' 35 The explanation of this interest- 
ing phenomenon is similar in all these situations. While 
for smooth particles (/? = —1) the rotational temperature 
is totally isolated from the translational one (so that the 
dotted arrows in Fig. [1] disappear), in the case of quasi- 
smooth particles (f3 > —1) a weak channel of energy 
transfer exists between the rotational and translational 
degrees of freedom and also the rotational temperature is 
subject to a weak cooling. Even though /3 might be very 
close to —1, the transfer of energy eventually becomes 



V. APPLICATION TO THE WHITE-NOISE 
THERMOSTAT 

As a second application, we consider now a homoge- 
neous granular gas subject to a stochastic thermostat 
force with properties of a Gaussian white noise £2r— The 
corresponding Boltzmann equation reads 4 ^ 



2 / Q 



where Xo i s a measure of the strength of the stochastic 
force. This force acts as a "thermostat" that injects en- 
ergy to the system, thus compensating for the collisional 
energy loss until a steady state is eventually reached. The 
evolution equations for the temperatures are 

d t T tr - mxo = -£ tr T tr , d t T Iot = _£ rot T rot , (5.2) 



(5.3) 



In terms of the dimensionless variables defined by Eqs. 
([3~T]) and ([53]) . Eq. (JTIJ) becomes 



a , . M20 - 3r d , , M02 9 

d '* + ^3— ^• (C0) + — ^- (W < 



r / d 



2 \dc 



J*[c,w|0], 



where 



(5.4) 



(5.5) 



-d s {(fw") +- [p ( M 20 - ST) + qm) (c p w q ) 

+ £p(p+l)( c P-2«,») = Mp9 . (5.6) 

In the steady state, Eq. (|5.2|) implies that £ tr T tr = mxo 
and £ rot = 0. Equivalently, 



ncr 2 (2T t Vm) 3 / 2 ' 
After taking moments in Eq. (|5.4[) one obtains 



M20 = 3r, 



M02 = 0. 



(5.7) 



(5.8) 
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FIG. 6: (Color online) Plot of the WNS temperature ratio 
T ro /T r vs the coefficient of normal restitution a for /3 = 
—0.9 (top panel), /3 = (middle panel), and /3 = 0.9 (bottom 
panel). The dashed and solid lines are the Maxwellian and 
Sonine approximations, respectively. 



As a consequence, Eq. (|5.6p yields, in the steady-state, 
5^20 = M40, (5.9) 



^20 — M22, 



/J04 = 0. 



(5.10) 



(5-11) 



Analogously to the HCS case, taking into account Eqs. 
(|3.18j) . (|3.19[) . and (|3.23p in the steady-state conditions 
(j5.8p - (|5.10p one obtains the cumulants 020, an, and ao2 
as functions of 9. Next, Eqs. (|3.25p and (|5.1ip provide a 
closed cubic equation for 9. 

Figure [5] shows the temperature ratio as a function of 
a for /3 = —0.9, 0, and 0.9. In the Maxwellian approx- 
imation [cf. Eq. p.lOp ] 9 is independent of a. On the 



0.06 



0.04 



0.02 



0.00 




0.06 
0.04 
0.02 
0.00 

0.02 
0.01 
0.00 



p=o 



p=-0.9 



P 0.9 

P- i ' 



P=-0.9 



P=0.9 



P=-l 



0.5 0.6 0.7 0.8 0.9 1.0 

a 

FIG. 7: (Color online) Plot of the WNS cumulants a 2 o (top 
panel), an (middle panel), and 002 (bottom panel) vs the 
coefficient of normal restitution a for /3 = — 1 (dotted lines), 
P = -0.9 (solid lines), f3 = (dashed lines), and /? = 0.9 
(dash-dotted lines). 



other hand, we observe that the Sonine approximation 
predicts a very weak dependence on a (note the verti- 
cal scales). Otherwise, the maximum deviation between 
both approximations in the range 0.5 < a < 1 is smaller 
than 1%, 2%, and 0.5% for /3 = -0.9, 0, and 0.9, re- 
spectively. It is clearly apparent that T rot /T tr increases 
with increasing roughness, ranging from in the smooth- 
sphere limit (/3 = —1) to about 1 in the opposite limit 
ft — > 1 . In the latter limit the Maxwellian approximation 
gives T rot /T tr = 1 for all a, while the Sonine approxima- 
tion gives T Iot /T tT ~ 1.01 for a = 0. 

The cumulants arc plotted in Fig.[7J Their magnitudes 
arc much smaller than in the HCS case (compare with 
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FIG. 8: (Color online) Plot of the WNS limit value 
lim^^_i(T rot /T tr )/(l + /3) vs the coefficient of normal resti- 
tution a. The dashed and solid lines are the Maxwellian and 
Sonine approximations, respectively. 



Fig. [3]). Apart from that, they are more significant for 
medium roughness (/? = 0) than for large (j3 = 1) or 
small (/? = —0.9) roughness. 

Again, it is worth analyzing separately the smooth- 
sphere limit j3 — >• — 1. From Eq. (|1.10[) we can expect 



(5.12) 



where •§ = 1 in the Maxwellian approximation. Thcrc- 
fore, in the limit (3 -> Eqs. (|3~Tg)) . 1(515)1 . and ([5^5)1 - 
(|3.25l) become 



1 + 
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(5.16) 



(5.17) 



Equations (|5.13|) and (|5 . 1 5[) include only the parameter 
a,2o an d so they are exactly the same as those obtained in 



the case of pure smooth spheres 4^ Therefore, application 
of the steady-state condition (|5.9p gives 



a-2o = 



-16(1 - 5) 
16(1 -a) 



1 - 85(1 - 5) 



239- 3275+ 1205 2 (2 - 5) 
1 -2a 2 



241 - 177a + 30a 2 (l-a)' 



(5.18) 



Thus, in contrast to the HCS case, the cumulant a 2 o in 
the WNS coincides in the limit j3 — s- — 1 with that at 
(3 = — 1. The Sonine approximation in that limit, Eq. 
(|5.18p , is also plotted in Fig. [7J where we observe that 
the curve is close to the one for f3 = —0.9, especially for 
large inelasticity. 

Application of Eqs. (|5"T3)) and (|5T6)) in Eq. (|5~TU)) im- 
plies an = 0. Next, from Eqs. ([5T5)| . ([Q3|) . and ([BTTgj) 
we get 



i) 



1 



3 n 

16 a 20 



1 — T§ a 2V) 

61 -45a + 6a 2 (l 



a) 



4[15-lla' + 2a 2 (l-a)]' 



(5.19) 



Finally, Eqs. ([5TT)) and ([5~T7| imply a 02 = 0. There- 
fore, in the limit (3 — > — 1 the translational distribution 
function is non-Maxwellian, as measured by a 2 o ^ 0, 
but the rotational distribution tends to a Maxwellian 
(an = ao 2 = 0) with a temperature T rot much smaller 
than T tr . Note that if the granular gas is made of strict 
smooth spheres (/3 = —1) the rotational distribution 
function (and hence the temperature T rot ) is not uniquely 
defined since it preserves its initial form. The parame- 
ter lim^_i(T rot /T tr )/(l +/3) = [k/2(1 + /t)]i? is plot- 
ted in Fig. [8) Again, the relative difference between the 
Maxwellian and the Sonine predictions is quite small (less 
than 0.7% in the range 0.5 < a < 1). 



VI. DISCUSSION AND CONCLUDING 
REMARKS 

The primary goal of this paper has been the derivation, 
within the Sonine approximation given by Eq. (|3.15l) . of 
the second- and fourth-degree collisional moments [cf. Eq. 
(|3.12j) ] in a granular gas made of inelastic rough hard 
spheres. The results are given by Eqs. (|3.18|) . (|3.19j) . and 
(|3.23[) - (|3.25p . In particular, the second-degree collisional 
moments (120 & n d fJ-02 are not but dimensionless versions 
of the collisional rates of change £ tr , £ rot , and £ associ- 
ated with the temperatures T tr , T Iot , and T, respectively 
[cf. Eqs. ([3"7T51) and (jSTTj)) ]. From that point of view, it 
is also worth emphasizing that Eqs. (|2.26)) (|2.28[) are ex- 
actly derived from the Boltzmann equation without any 
further assumption, so that they are not restricted to any 
Sonine approximation. 

Our results represent extensions of some previously de- 
rived results. On the one hand, Eqs. (|3.18)) - (|3.20j) are So- 
nine extensions of those obtained in the Maxwellian ap- 
proximation defined by Eq. (|1.2[) j 29 i 34 On the other hand, 
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Eqs. (|3.18[) and (|3.23|) are extensions to rough spheres of 
previous Sonine derivations for smooth spheres^ 

Since the Boltzmann collision operator [cf. Eq. (|2.20l) ] 
is local in time and space, Eqs. p.!8[) - p.25[) keep being 
applicable to inhomogeneous and unsteady states. They 
also hold in the context of the Enskog equation for homo- 
geneous states, except that the collision frequency (|1.7[) 
must be multiplied by the pair correlation function at 
contact, g(a). Apart from that, it is important to bear 
in mind that some restrictions apply to the Sonine ap- 
proximation p,15|) . First, it has been assumed that the 
mean angular velocity vanishes, i.e., (u>) = 0. This re- 
striction is easy to circumvent by replacing u) — > u> — (u>) 

and T rot T™* [see discussion below Eq. (f2~25"l) ] in Eqs. 
Q, (|3~Tj) . and (|3~2"j) . These changes would affect the 
collision rules in Eqs. p.3p - (|3.5p by the changes 9 — > 

T rot /T tr and w x + w 2 -»■ Wl + w 2 + 2(w)/(2T rot //) 1/2 , 
given that the angular velocity is not a a conserved quan- 
tity. For the expressions of £ tr and £ rot in the Maxwellian 
approximation with (lj) =/= 0, the reader is referred to 
Refs.HH 

As a second restriction, notice that Eq. p,15[) may 
be less useful in strongly anisotropic states where the 
dependence of 0(c, w) on the six velocity components is 
not exhausted by the three scalar quantities c 2 , w 2 , and 
(c • w) 2 . Finally, while Eq. p,15p treats the two second- 
degree moments as independent quantities, it does not 
do so with the a priori four independent fourth-degree 
moments. Instead, Eq. (|3.15p assumes that the moment 
((c • w) 2 ) is enslaved to (c 2 w 2 ) by Eq. p,17p . Therefore, 
the study of the orientational correlation between the 
translational and angular velocities has not been included 
in our scheme. 



We have applied Eqs. p.!8p - p.25p to two paradig- 
matic homogeneous and isotropic situations: the sim- 
ilarity solution of the HCS and the steady-state solu- 
tion of the WNS. In both cases we have found that 
the Maxwellian approximation for the temperature ratio 
T rot /T tr , being much simpler than the corresponding So- 
nine approximation, does a very good job, especially for 
large or small roughness. On the other hand, departures 
of the velocity distribution function from the Maxwellian, 
as measured by the cumulants 020, an, and 002, cannot 
be ignored. This is especially true in the case of the HCS 
for medium and small roughness. In fact, in the quasi- 
smooth limit P — > — 1 the HCS results differ markedly 
from those obtained in the case of pure smooth spheres 
(/3 = —1). This interesting singular behavior is directly 
related to the unsteady character of the HCS and thus it 
is absent in the steady WNS. 

We expect that this work can contribute to our under- 
standing of the subtle interplay between roughness and 
inelasticity in granular gases and how the former feature 
modifies the properties of inelastic smooth spheres. We 
plan to assess the qualitative and quantitative results de- 
rived here by comparison with computer simulations for 
several situations of physical interest. 

Acknowledgments 

This paper is dedicated to the memory of Carlo Ccrcig- 
nani. The work of A.S. has been supported by the Minis- 
terio de Ciencia e Innovacion (Spain) through Grant No. 
FIS2010-16587 (partially financed by FEDER funds). 
The work of G.M.K. and M.S. has been supported by 
the Conselho Nacional de Desenvolvimento Cienti'fico c 
Tecnologico (Brazil). 



* Electronic address: [andres@unex.es 

URL: http : //www.unex. es/eweb/f isteor/andres/ 
t Electronic address: kremer@fisica.ufpr.br 

1 J. A. Carrillo, C. Cercignani, and I. M. Gamba, "Steady 
states of a Boltzmann equation for driven granular media," 
Phys. Rev. E 62 7700 (2000). 

2 C. Cercignani, "Shear flow of a granular material," J. Stat. 
Phys. 102, 1407 (2001). 

3 C. Cercignani, R. Illner, and C. Stoica, "On diffusive equi- 
libria in generalized kinetic theory," J. Stat. Phys. 105, 
337 (2001). 

4 C. Cercignani, "The Boltzmann equation approach to the 
shear flow of a granular material," Phil. Trans. R. Soc. 
Lond. A 360, 407 (2002). 

5 A. V. Bobylev and C. Cercignani, "Moment equations for 
a granular material in a thermal bath," J. Stat. Phys. 106, 
547 (2002). 

6 A. V. Bobylev and C. Cercignani, "Self-similar asymptotics 
for the Boltzmann equation with inelastic and elastic in- 
teractions," J. Stat. Phys. 110, 333 (2003). 

' A. V. Bobylev, C. Cercignani, and G. Toscani, "Proof of an 
asymptotic property of self-similar solutions of the Boltz- 
mann equation for granular materials," J. Stat. Phys. Ill, 



403 (2003). 

8 A. V. Bobylev, C. Cercignani, and I. M. Gamba, "Gener- 
alized Kinetic Maxwell Type Models of Granular Gases," 
in Mathematical Models of Granular Matter, edited by G. 
Capriz, P. Giovine, and P. M. Mariano, Lecture Notes in 
Mathematics, Vol. 1937 (Springer, Berlin, 2008), pp. 23-57 

9 A. V. Bobylev, C. Cercignani, and I. M. Gamba, "On 
the Self-Similar Asymptotics for Generalized Nonlinear Ki- 
netic Maxwell Models," Comm. Math. Phys. 291, 594 
(2009). 

10 I. Goldhirsch, "Rapid granular flows," Annu. Rev. Fluid 
Mech. 35, 267 (2003). 

11 S. J. Moon, J. B. Swift, and H. L. Swinney, "Steady-state 
velocity distributions of an oscillates granular gas." Phys. 
Rev. E 69, 011301 (2004). 

12 J. T. Jenkins and M. W. Richman, "Kinetic theory for 
plane flows of a dense gas of identical, rough, inelastic, 
circular disks," Phys. Fluids 28, 3485 (1985). 

13 C. K. K. Lun and S. B. Savage, "A Simple Kinetic Theory 
for Granular Flow of Rough, Inelastic, Spherical Particles," 
J. Appl. Mech. 54, 47 (1987). 

14 C. S. Campbell, "The stress tensor for simple shear flows 
of a granular material," J. Fluid Mech. 203, 449 (1989). 



14 



15 C. K. K. Lun, "Kinetic theory for granular flow of dense, 
slightly inelastic, slightly rough spheres," J. Fluid Mech. 
233, 539 (1991). 

16 C. K. K. Lun and A. A. Bent, "Numerical simulation of 
inelastic spheres in simple shear flow," J. Fluid Mech. 258, 
335 (1994). 

17 C. K. K. Lun, "Numerical simulation of inelastic spheres 
in simple shear flow," Phys. Fluids 8, 2868 (1996). 

18 P. Zamankhan, H. V. Tafreshi, W. Polashenski, P. Sarko- 
maa, and C. L. Hyndman, "Shear induced diffusive mix- 
ing in simulations of dense Couette flow of rough, inelastic 
hard spheres," J. Chem. Phys. 109, 4487 (1998). 

19 J. T. Jenkins and C. Zhang, "Kinetic theory for identical, 
frictional, nearly elastic spheres," Phys. Fluids 14, 1228 
(2002). 

20 W. Polashenski, P. Zamankhan, S. Makiharju, and P. 
Zamankhan, "Fine structures in sheared granular flows," 
Phys. Rev. E 66, 021303 (2002). 

21 I. Goldhirsch, S. H. Noskowicz, and O. Bar-Lev, "Nearly 
Smooth Granular Gases," Phys. Rev. Lett. 95, 068002 
(2005). 

22 A. Goldshtein and M. Shapiro, "Mechanics of collisional 
motion of granular materials. Part 1. General hydrody- 
namic equations," J. Fluid Mech. 282, 75 (1995). 

23 M. Huthmann, and A. Zippelius, "Dynamics of inelasti- 
cally colliding rough spheres: Relaxation of translational 
and rotational energy," Phys. Rev. E 56, R6275 (1998). 

24 S. McNamara and S. Luding, "Energy non-equipartition 
in systems of inelastic, rough spheres," Phys. Rev. E 58, 
2247 (1998). 

25 S. Luding, M. Huthmann, S. McNamara, and A. Zippelius, 
"Homogeneous cooling of rough, dissipative particles: The- 
ory and simulations," Phys. Rev. E 58, 3416 (1998). 

26 O. Herbst, M. Huthmann, and A. Zippelius, "Dynamics 
of inelastically colliding spheres with Coulomb friction: 
Relaxation of translational and rotational energy," Gran. 
Matt. 2, 211 (2000). 

27 T. Aspelmeier, M. Huthmann, and A. Zippelius, "Free 
Cooling of Particles with Rotational Degrees of Freedom," 
in Granular Gases, edited by T. Poschel and S. Luding 
(Springer, Berlin, 2001), pp. 31-58. 

28 R. Cafiero, S. Luding, and H. J. Herrmann, "Rotationally 
driven gas of inelastic rough spheres," Europhys. Lett. 60, 
854 (2002). 

29 A. Zippelius, "Granular gases," Physica A 369, 143 (2006). 

30 S. Luding, Phys. Rev. E "Granular materials under vibra- 
tion: Simulations of rotating spheres," 52, 4442 (1995). 

31 N. Mitarai, H. Hayakawa, and H. Nakanishi, "Collisional 



Granular Flow as a Micropolar Fluid," Phys. Rev. Lett. 
88, 174301 (2002). 

32 N. V. Brilliantov, T. Poschel, W. T. Kranz, and A. Zip- 
pelius, "Translations and Rotations Are Correlated in 
Granular Gases," Phys. Rev. Lett. 98, 128001 (2007). 

33 W. T. Kranz, N. V. Brilliantov, T. Poschel, and A. Zip- 
pelius, "Correlation of spin and velocity in the homoge- 
neous cooling state of a granular gas of rough particles," 
Eur. Phys. J. Spec. Top. 179, 91 (2009). 

34 A. Santos, G. M. Kremer, and V. Garzo, "Energy produc- 
tion rates in fluid mixtures of inelastic rough hard spheres," 
Prog. Theor. Phys. Suppl. 184, 31 (2010). 

35 A. Santos, "Homogeneous Free Cooling State in Binary 
Granular Fluids of Inelastic Rough Hard Spheres," in Rar- 
efied Gas Dynamics: Proceedings of the 27th Interna- 
tional Symposium on Rarefied Gas Dynamics, D. A. Levin, 
ed. (AIP Conf erence Proceeding s, Melville, NY, 2011), in 
press; preprint arXiv:1007.0701 

36 A. Santos, "A Bhatnagar-Gross-Krook-like Model Kinetic 
Equation for a Granular Gas of Inelastic Rough Hard 
Spheres," in Rarefied Gas Dynamics: Proceedings of the 
27th International Symposium on Rarefied Gas Dynamics, 
D. A. Levin, ed. (AIP Conference Proceedings, Melville, 
NY, 2011), in press; preprint [arXiv:1007.0700l 

37 H. Uecker, W. T. Kranz, T. Aspelmeier, and A. Zippelius, 
"Partitioning of energy in highly polydisperse granular 
gases," Phys. Rev. E 80, 041303 (2009). 

38 D. R. M. Williams and F. C. MacKintosh, "Driven gran- 
ular media in one dimension: correlations and equation of 
state," Phys. Rev. E 54, R9 (1996). 

39 D. R. M. Williams, "Driven granular media and dissipa- 
tive gases: correlations and liquid-gas phase transitions," 
Physica A 233, 718 (1996). 

40 M. R. Swift, M. Boamfa, S. J. Cornell, and A. Maritan, 
"Scale invariant correlations in a driven dissipative gas," 
Phys. Rev. Lett. 80, 4410 (1998). 

41 T. P. C. van Noije and M. H. Ernst, "Velocity distributions 
in homogeneous granular fluids: the free and the heated 
case," Gran. Matt. 1, 57 (1998). 

42 J. M. Montanero and A. Santos, "Computer simulation 
of uniformly heated granular fluids," Gran. Matt. 2, 53 
(2000). 

43 A. Santos and J. M. Montanero, "The second and third So- 
nine coefficients of a freely cooling granular gas revisited," 
Gran. Matt. 11, 157 (2009). 

44 J. J. Brey, J. W. Dufty, and A. Santos, "Dissipative dy- 
namics for hard spheres," J. Stat. Phys. 87, 1051 (1997). 



